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We analyze pairing in two dimensional spin liquids. We argue that interplane pairing enhanced 
by magnetic correlations is the most plausible explanation of the spin gap phenomenon observed in 
underdoped cuprates. The details of the pairing theory depend on the in-plane antiferromagnetic 
correlations. We consider two models: 2pF correlations induced by a strong gauge field interaction 
and undamped spin waves. We estimate the pairing temperature T s and the angular dependence of 
the gap function and discuss physical consequences. 



I. INTRODUCTION 

Many high T c materials exhibit anomalous tempera- 
ture dependence of the bulk magnetic susceptibility, \i 
in a range of temperatures above T c , suggesting that a 
spin pseudogap opens above the superconducting tran- 
sition temperaturdl. To define more precisely what we 
mean by spin pseudogap consider the susceptibility data 
for Y Bd2Cu^O% shown as the curve labeled "248" in Fig. 
1. (The data points shown in Fig. L-jwere obtained 
from tlie susceptibility reported in RefEfl as discussed 
in Ref.cl) There are clearly two regimes, separated by a 
scale T s « 200K. For T > T s , X = A + BT with A, 
B > 0. For T < T s , x{T) drops more rapidly; indeed a 
straight line fit to X (T) for T c < T < T s (T c = 80K is 
the superconducting transition temperature for this com- 
pound) would yield a negative x a t T = 0. A negative 
x(0) is impossible; therefore the negative extrapolation 
implies that even for T > T c (but T <T S ) there is a gap 
for spin excitations. The origin of this "spin" gap is the 
main focus of this paper. 

The properties of high temperature superconducting 
materials are anomalous, and many different theories 
have been proposed to describe them. However, no gen- 
eral consensus has emerged. The_j^igia. of. spin gap 
has been discussed by many authorsBd^Qu'LTU but these 
treatments are not completely satisfactory because, as we 
shall argue below, they are based on models which do not 
agree with all available data. We believe that any theory 
of the spin gap should have the following ingredients: 

(i) At least some of the magnetic response is Pauli-like, 
i.e. it comes from a particle-hole continuum of spin 
1/2 fermions. 

(ii) Formation of the spin gap involves pairing instabil- 
ity of these fermions. 



(iii) This pairing is not a superconducting pairing, i.e. 
does not produce a Meissner effect or paraconduc- 
tivity. 



These assumptions imply that paired fermions are neu- 
tral and that any theory of the spin gap must involve the 
phenomenon of spin-charge separation; we use the gauge 
theory formalism to describe this. 
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FIG 1 . Susceptibility of several high T c materials: one- 
plane La,2- x Sr x CuO<i with varying doping x which, as 
we shall argue, does not show spin gap behavior and bi- 
layer Y BaiCu/^D% (248) which shows it. Susceptibility 
is measured in the units of fj, 2 B eV per Cu atom. 



Points (i)-(iii) do not completely specify the model. 
There are two additional issues. 
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(iv) It is widely believed that there are strong antiferro- 
magnetic spin fluctuations in high-T c materials. In 
the body of the paper we summarize the evidence 
for the existence of antiferromagnetic spin fluctu- 
ations, outline the different theories proposed to 
describe them, and give the implications for our 
calculations. 

(v) We believe that bilayer structure of Y Ba^Cu^O^ 
is important for the formation of the spin gap be- 
cause it leads to interplane pairing. Whether the 
spin gap exists in single-plane materials such as 
La2~ x Sr x CuOi, in which interplane coupling is 
geometrically frustrated, is controversially. Cer- 
tainly, the evidence for it is weaker (see Fig. 1). 
Our theoretical results imply that spin gap effects 
are greatly enhanced in bilayer or multilayer sys- 
tems. 

The outline of the paper is as follows. In Section II, we 
review the relevant experimental data and show that it 
implies points (i)-(v) above. In Section III we formulate 
the theoretical model and identify two different scenar- 
ios. In Section IV and V we present the solutions for the 
different scenarios. Section V contains a comparison of 
results to data and a conclusion. 



II. DATA 

(i) Particle-hole continuum: We begin with the evidence 
that at least some of the magnetic response in spin-gap 
systems such as Y Bo2CuaO^ comes from a particle- hole 
continuum of spin 1/2 fermions. A phenomenological ar- 
gument is that there is little doubt that such a contin- 
uum describes the magnetic properties of YBa 2 Cu^O-j. 
The uniform susceptibility of this compound has magni- 
tude slightly larger than predicted by band theory and 
has very weak temperature dependence (as expected for 
a Fermi liquid)!! 2 ]. The oxygen NMR relaxation rate 
\/ 17 TiT is temperature independent (as expected in a 
Fermi liquid) and has magnitude, slightly larger than pre- 
dicted by the Korringa relationll3. The copper relaxation 
rate l/ 63 TiT increases as the temperature is decreasedliJ; 
this has been argued to be due to antiferromagnetic cor- 
relations within a Fermi liquid statec3. Also, photoe- 
mission experimental^ have observed a large Fermi sur- 
face in optimally doped Y Bo2Cu^Ot. The spin-gap com- 
pounds YBa2CusOe.e and Y Ba2Cu^O^ are produced by 
removing carriers from Y Ba2Cuj > 0-j\ further, the physi- 
cal properties seem to vary smoothly with doping. There- 
fore, it is natural to assume that the particle hole con- 
tinuum is also present in the spin-gap compounds. 

An alternative point of viewu is that the important 
spin excitations in underdoped compounds are weakly 
damped, gapped antiferromagnetic spin waves, with dis- 
persion lo 2 = c 2 k 2 + A 2 . However, theories in which 
the only magnetic excitations are spin waves with a gap 



are inconsistent with the Xs(T) data at high tempera- 
tures (T > T s ), because at temperatures greater than 
the value of the zero-temperature gap A , these theories 
predictO \s = A+BT with A < if A > and also pre- 
dict for YBa2Cu40*-& value for B larger, by a factor of 
six, than is observecQ Another difficulty with the spin- 
wave-only model is that the kinematics of spin waves im- 
plies that the oxygen relaxation rate drops very rapidly 
with temperature. Indeed the theoretical prediction is 
1/piTxs] ~ T in contrast to the experimental result 
l/(TiTxs) ~ constant$3 Therefore there is little doubt 
that one must at least supplement the spin wave model 
with a particle-hole continuum ofJsawiions which controls 
at least the small q spin responsdl3E3. We note, however, 
that although the existence of such a particle-hole contin- 
uum would increase the constant part, A, of the uniform 
susceptibility and would add a constant contribution to 
the oxygen relaxation rate, it would not affect the slope, 
B, of the susceptibility if the interaction between spin 
waves and particle-hole continuum is weak. There is as 
yet no model which describes all of the data. In any 
event, to understand the spin gap observed in Xs{T) one 
must understand how to open a gap in the spectrum of 
the fermions making up the particle-hole continuum. 

(ii) Singlet pairing: In the absence of charge density wave 
order the only possibility for suppressing the spin re- 
sponse of fermions is to pair them into singlets. One ex- 
ample of such pairing is the BCS superconducting state. 
The pairing scale T s 200 K is much less than any mi- 
croscopic scale such as J or the bandwidth, so the pairing 
must be understood as a low energy instability of the sys- 
tem. 

(iii) Superconductivity: The data strongly suggest that 
the spin gap is not due to incipient superconducting pair- 
ing. There are two arguments; one is that the spin gap 
scale T s in Y Ba2Cu^O% is approximately 200K, far above 
the maximum T c ss 95K observed in any member of the 
YBCO family. The second is that superconducting fluc- 
tuations normally have a dramatic effect on the resis- 
tivity, which begins to drop as T is decreased through 
the scale at which the fluctuations begin and drops more 
and more rapidly as T c is approached. This behavior is 
not observed in Y Ba,2CuzO§.§ and in Y Ba 2 Cu±0%. In 
these materials, there is some drop in p as T is decreased 
through T s , but then a flattens out and depends only 
weakly on T near 7 1 J 21 tr 2 l The drop in p for T ~ T s 
has been attributecBa, in our view correctly, to changes 
in the inelastic scattering mechanism associated with the 
onset of the spin gap. It has been argued that negative 
U Hubbard model, in which the important physics is sin- 
glet pairing of conventional electrons describes the spin 
gap phenomenaO. In our view the strong evidence for re- 
pulsive interactions and against paraconductivity renders 
this model irrelevant. 

To summarize: the data imply the existence of a 
particle-hole continuum of spin excitations and require 
that these excitations be paired into singlets in a way 
that does not produce superconductivity. This implies 
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that "spin-charge separation" must occur. There are 
many scenarios of spin-charge separation, all stemming 



from Anderson's origi 
gauge theory approacl 



topcEa^. We shall adopt the 
We also note that unlike the 
conventional superconducting pairing of electrons, the 
pairing of chargcless fermions does not necessary imply 
the breaking of any symmetry, and so max result in a 
crossover, not in a genuine phase transitionB. 
(iv) Antiferromagnetism: A large literature has devel- 
oped around the issue of antiferromagnetic correlations 
in high T c materials. Many controversies remain unre- 
solved, but it is generally agreed that NMR T\ and Ti 
experiments imply the existence of strong, temperature 
dependent antiferromagnetic correlations. Specifically, 
for YBCI2CU4O8, the two Cu relaxation rates, 1/T\T 
and 1 /T^j increase roughly as 1/T as T is decreased (for 
T > T s )E3'tZl. Generally, the NMR rates are given by 
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x"M 



TiT 
1 

% 



1 1/2 



(1) 
(2) 



where F q is a form factor which is different for different 
nuclei. 

The only tenable interpretation of the temperature 
dependence of the copper relaxation rate is that both 
real and imaginary parts of x( ( l) diverge at a particular 
wavevector Q, i.e. that there is an incipient magnetic 
instability. To prove this, suppose on the contrary that 
the T-dependence of l/(TiT) came from a wide range of 
momenta, so cj) = tfi(q)f(u>,T), with 4>(q) a tem- 
perature independent function. For large frequencies, 
lu ^> T, f(oj) should not depend on temperature, so in 
particular the imaginary part f"(ui) — Aiu x , whereas 
at small frequencies, u> <C T, f"(to) should be propor- 
tional to frequency, f"(uj) — A2L0. The proportional- 
ity coefficient, A2, can be estimated by matching the 
low frequency formula and the high frequency formula 
at uj ss T, yielding A 2 = A X T X . The NMR T x data 
imply x ~ which via the Kramers-Kronig relation im- 
plies x'(<?,cj = 0) ~ InT. This temperature dependence 
is too weak to account for the T2 data, so the hypoth- 
esis of a momentum-independent divergence of x must 
be rejected. It should be noted, however, that although 
neutron scattering experiments detect antiferromagnetic 
fluctuations, neutnon and NMR data are not quantita- 
tively consistentE§E2l 

The proper theoretical model for the antiferromagnetic 
fluctuations is not clear. The two principal proposals 
are that the dominant antifaeromagnetic excitations are 
weakly damped spin wavesaO, or are particle-hole pairs 
of an antiferromagnetically correlated fermion systemLj. 
The weakly damped spin wave picture applies to the 
magnetic insulating parent compound and by continuity 
might be expected to apply to lightly doped but non- 
ordered materials. The particle-hole picture presumably 



applies to the optimally doped materials (which have 
been shown by photoemission to have large (Luttinger) 
Fermi surface) and by continuity might be expected to 
apply to somewhat underdoped materials. The crossover 
between these two regimes is a active area of researched 
but has not been understood in detail. We consider im- 
plications of both pictures for the pairing interaction. 

The model of particle-hole pairs requires further dis- 
cussion. We argued above on the basis of resistivity data 
that the fermions must be charge zero objects. An addi- 
tional argument against conventional charge e Fermi liq- 
uid with antiferromagnetic correlations is that the mag- 
netic properties of conventional ferrninUauid, are incom- 
patible with the Ti, T 2 and X s dataBEZHM The only 
known model of spin 1/2 charge fermionic excitations 
("spinons") which does not break time reversal or par- 
ity symmetry is the "spin liquid" . In a spin liquid, the 
gauge interaction between fermions affects the magnetic 
properties and has been shown to lead to a reasonable 
description of the spin dynamics of high T c materialslEI 
In the spin-wave picture a particle-hole continuum is also 
present; for the reasons given above we must assume that 
the underlying fermionic excitations are spinons. How- 
ever, in this case the spin waves dominate the large q 
magnetic response. 

(v) Bilayers: Our discussion so far has been focussed on 
members of the Y — Ba family of high T c materials. We 
believe that the magnetic dynamics of members of other 
families of high T c materials are similar in all respects ex- 
cept for the existence of a spin gap in a wide temperature 
range otf^e T c . It is possible that, as some authors have 
arguedtrLij, La2- x Sr x CuO '4 shows the beginning of spin 
gap behavior for T ~ 60K, relatively near the supercon- 
ducting T c . However, to our knowledge only underdoped 
members of the Y — Ba family exhibit spin gap behavior 
over a wide range of temperatures above T c . 

Consider first the \s f° r La2- x Sr x CuC>4 data pre- 
sented in Fig. 1. These are obtained from bulk suscepti- 
bility data biz-subtracting core and van-Vleck susceptibil- 
ities given inE.1 Note that although there is a downturn in 
Xs(T) at T < 100K, extrapolated value of limx->o Xs(T) 
is positive and relatively large. Although uncertainties in 
the value of the van-Vleck susceptibility may exist, these 
are in our view by no means laxge enough tq-praduce a 
limr->oXs(T) < 0. The oxygenEfl and coppert3'E£l T\ re- 
laxation data similarly show no sign of an extra downturn 
at a temperature T > T c (although the rate of increase 
of the copper Tj -1 may slow for T < 80 K). Note that at 
T > T s « T c , La 2 ~ x Sr x Cu04 does exhibit \ s ss A + BT 
regime as well as other properties difficult to describe 
in either a Fermi liquid or a purely spin wave picture. 
We infer from this that the bilayer structure of the Y-Ba 
material is important only for raising T s sufficiently far 
above T c that spin gap effects are easily observable. 

There is in fact substantial evidence that the spin 
degrees of freedom in different planes of a bilayer 
are strongly coupled. Neutron scattering measure- 
ments have essentially only detected spin fluctuations 
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in which spin_.oa .adjacent Cu0 2 layers are perfectly 
anticorrelatec£3t3~E3. Moreover, the coupling between 
Cu spins on adjacent planes has been directly measured 
in a recent NMR T 2 experiment, the Cu nuclear spins 
in one plane of a bilayer were pumped and Cu nuclear 
spins in the other were measurecO. This experiment de- 
termines the cross relaxation time T 2 * which is given by 
expressiorO 



£(^x'i 2 M)) 5 



1/2 



(3) 



This expression is very similar to the expression (Q) for 
in-plane relaxation rate, except that instead of a single 
plane susceptibility x'(g, 0) it contains Xizili^) which 
measures the response of spins on plane 1 to the mag- 
netic field on plane 2. 

We assume that electrons on adjacent planes in bilayer 
interact via the Hamiltonian 



(1)5(2) 



(4) 



If the interaction, Jj_ is weak then the between the planes 
susceptibility, X12, m &y be calculated by perturbation 
theorycil and is 



xi 2 (g.o) = J±[x'(q,0)f 



(5) 



Experimentally the ratio T 2 /T 2 * grows from .14 at 200K 
to 0.28 at 100K. This increase reflects the temperature 
dependence of \' ■ If x' is divergent at some wavevec- 
tor Q as T — > and is given by a scaling form x' = 
T~ a f{\q — Q\T~ X ) (where x and a are scaling expo- 
nents) then from Eqs. (|]||,|^) it may be shown that 
T 2 /T 2 * — Cr-J±x'{0,Q) with c a constant of the or- 
der of unityHil. Thus, the observed maximal ratio of 
T2/T2 — 0.3 implies that the interplane coupling is not 
negligible, but still may be treated via perturbation the- 
ory. In this paper we shall show that the effect of this 
interplane coupling on the fermions is large and in fact 
leads to the opening of a spin gap. 

Eq. (Q) applies to YBaCuO and to the multilayer 
BSCCO compounds in which the Cu ion in one plane 
sits directly over the Cu ion in the next lower plane. It 
does not apply to La 2 - x Sr x CuO^ compounds in which 
the crystal structure is such that a Cu ion on plane is 
coupled equally to four Cu in each adjacent plane, so 
Eq. (0) would become 



E,i q(a) q(<>.+ l) 
J ij°i °; 



(0) 



The crystal structure of La 2 - x Sr x CuO^ implies that 
J±(q) vanishes at q — (71", 7r), so the enhancement of inter- 
plane pairing by antiferromagnetic fluctuations is much 
less effective than in the YBaCuO or BSCCO systems. 



III. MODEL 

(i) Single plane. In Section II we showed that exper- 
iment implies that a theoretical treatment of spin gap 
effects in high-T c superconductors should involve pairing 
of fermions in a spin liquid. In this subsection we de- 
scribe the model we use for the spin liquid in one CUO2 
plane liquid and discuss different pairing mechanisms. 

The low energy excitations of a spin liquid are S = 1/2, 
charge fermions, c 1 , near a Fermi line, and a bosonic 
gauge field, a. The action describing the spin liquid has 
been derived from more fundamental models of correlated 
electrons by assuming that spin charge separation exists, 
i.e. that the electron field ip + may be written as the 
product of a spinless, charge e Bose field b and a, S = 1/2 
charge Fermi field c + as = c + b and that the effect of 
the charge degrees of freedom on the spin degrees of free- 
dom is small. These assumptions have been shown to be 
justified in the low -doping, large spin degeneracy limit 
of the t — J modelcfl. Whether these assumptions are 
theoretically justifiable in the physically relevant regime 
is still controversial. We shall assume that they are be- 
cause we see no other way to explain the experimental 
data discussed in Section II. 

The spin liquid model is specified by the action 

+£ Wcl ua Cp 2 , a cl St „Cp it<T 8(£2pi) + -T-nfl u (7) 
Pi 9a 

Here, as usual, = d^a v — d^a^, go is the bare fermion- 
gauge field interaction constant, v = W is a constant 
of short-range interaction and a = 1 ... N is a spin in- 
dex. In the physical case the spin degeneracy N = 2. 
The gauge field modifies the properties of the fermions. 
These modifications have been studied in detail. Thp r.e-. 
suits we shall need are: (i) the electron self energy ial3c3 
S(e) = uj^ 3 e 2 ^ 3 (ii) the vertex, T, coupling the fermion 
spin to an external magnetic field of wavevector q be- 
comes singular at |g| = 2pp while at allxtther wavevectors 
the vertex corrections are not singulars. Specifically, 



(8) 



Here A is the largest of «f(M - 2p F ), uj 2 / 3 , T 2 / 3 and 



w 



2V3 



334 
v f9o 

ir 2 p Q 



(9) 



is the upper cut-off scale determined by the strength 
of the gauge field fluctuations. By "2pp" we mean a 
wavevector Q which connects two points on the Fermi 
line with parallel tangents. For a circular Fermi line any 
vector Q of magnitude 2p F connects two such points. 
The exponent a has been calculated only in the limits 
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N ^> 1 and N <C 1 , by extrapolation, these results to the 
physical value N = 2 we estimatecQ that 1 > a > 1/3. 

The main effect of the self-energy renormalization is 
that the resulting inelastic scattering rate ~ T 2 / 3 is so 
strong that no non^singular interaction can lead to a BCS 
pairing of spinonaij (except via a first order transition 
which is not observed); therefore any spinon-based the- 
ory of the spin gap must involve a singular interaction. 

Two cases arise for the vertex renormalization. If 
a < 1/3, the spin physics is not modified in an essential 
way. A T = critical point separates a W < W c phase 
with short range spin correlations from a W > W c phase 
with long range order; the appearance of the anomalous 
exponent a in the "2pp" vertex modifies the critical prop- 
erties of-,the transition at W = W c as discussed in detail 
in RefM However, if a > 1/3 then the < W < W c 
phase is anomalous, and exhibits a divergent ll 2pp" spin 
susceptibility and power law correlations 



X(^,k) 



uoPo 



1 





f\u\\ 


Cuj 





2(7-2/3 



/ |fc|||uF \ 



3cr-l 



(10) 



where c w ,Cfc ~ 1. Here we use local momentum co- 
ordinates associated with the Fermi line, namely k = 
Q + en fcii where Q connects two points on the Fermi line 
with parallel tangents and eii is the unit vector parallel 
to the Fermi velocity at these points. These coordinates 
are generalization of radial fen = \p\ — 2pp and angular co- 
ordinates for the case of non-circular Fermi line. The ex- 
perimental implications have been discussed elsewhereo; 
in particular, it has been shown that the choice a = 3/4 
yields rough agreement with experimental data. 

Having discussed the spin liquid model we now con- 
sider possible pairing interactions. The short range in- 
teraction W would lead to pairing in a Fermi liquid; how- 
ever, as mentioned above, for a spin liquid the inelastic 
scattering due to the gauge field suppresses any second- 
order pairing instability due to a non-singular interac- 
tion!-. ^ or ^ n ' s reason we believe the results obtained in 
Ref.El do not explain the spin gap phenomena— Of course, 
a first order transition would still be possibles; but there 
is no experimental evidence for this in high T c materials. 
However, singular interactions exist. One involves the 
gauge field, but this interaction is repulsive in all chan- 
nels for— the model specified above and does not lead to 
pairingcfl Another singular interaction comes from ex- 
change of long-ranged spin fluctuations; these may arise 
either from proximity to a T — antiferromagnetic tran- 
sition or because a > 1/3. 

There are three types of antiferromagnetic transitions, 
distinguished by the relation of the ordering wavevector, 
Q, to 2pp. If |Q| > 2pp, the fermion-fermion interac- 
tion mediated by spin fluctuations is in fact not singular 
for the fermions near the Fermi line. If |Q| = 2pp and 
a < 1/3, then the singularity of the interaction is too 
weak to overcome the pairbreaking effect of the gauge 



field. If Q < 2pp, then one obtains a logarithmic diver- 
gence in the pairing amplitude. The theory of the pairing 
for Q < 2pp case mayJae derived by the following the ar- 
guments given in Ref.EUl but replacing the factor x 2 ( w i l) 
by the first power x(k>, qjUn Eq. (5) of Ref.E2l. The steps 
leading to Eq. (8) of Ref.EHl yield a logarithmic divergence 
of the pairing kernel in the gap equation. However, we do 
not believe that the Q < 2pp case is relevant to high T c 
materials because the predicted temperature dependence 
of the copper NMR T 2 rate is too weak and because neu- 
tron scattering has oniw-observed fluctuations peaked at 
wavevectors Q > 2pp 

If a > 1/3, then the susceptibility is divergent 
Q — ► 2pp, (uj,T) — ► 0. However, as was shown in Re 
the same physics implies that the fermion-fermion inter- 
action is renormalized to zero (due to the renormaliza- 
tion in the Cooper channel), so the divergence in the 
susceptibility does not propagate into any other physical 
quantity. 

(ii) Two Planes: The theoretical discussion in the previ- 
ous subsection and the experimental analysis of Section II 
implied that theories involving only a single CuOi plane 
could not explain the existence of a spin gap in a wide 
range of temperatures above T c . In this subsection we 
extend the theory of the spin liquid to include interplane 
coupling. We assume that each plane is described by 
the Hamiltonian, Eq. (R), and that the only coupling be- 
tween the planes is given by Eq. (^). In particular, terms 
of the form t^c + ^c^> + h.c. are not allowed: there is no 
coherent hopping of spinons between planes. The as- 
sumption that there is no interplane hopprng of spinons 
may be justified by extending the derivationEj of the spin- 
liquid action, H, to the two plane situation. In the micro- 
scopic model there is a between-planes electron hopping 
t e l_ which, it is reasonable to assume, is much less than 
the in-plane hopping, t. It emerges from the theoretical 
derivation that one must have t e ^_ greater than a critical 
value of order t in order to have coherent between-planes 
hopping of spinons. If there is no coherent hopping then 
the leading coupling term is Eq. (^). 

We assume that J± is sufficiently weak that it may be 
treated in perturbation theory. This assumption is justi- 
fied by the cross-relaxation experiment of Stern et aH, 
as previously discussed. The only effect wc need con- 
sider the pairing interaction due farJ-f- If Q < 2pp then 
the relevant theory was given inliixEI If Q = 2pp and 
<7 < 1/3 the theory is very similar. As discussed above, 
we do not believe any of these starting points are con- 
sistent with experiment. We therefore study in Section 
IV the case a > 1/3. In Section V we treat the case of 
undamped spin fluctuations. 



IV. GAP EQUATION; a > 1/3 



In this section we consider pairing due to the between 
planes interaction, Eq. (ji|). In general, a pairing insta- 
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bility of a fermion system is signalled by a divergence of 
the series .of particle-particle ladder diagrams as shown 
in Fig. 2.t3 In the present problem, the two lines cor- 
respond to fermions on different planes and the between 
planes interaction J± is renormalized by the gauge inter- 
action. This rcnormalization implies that the basic rung 
of the ladder is an effective pairing interaction 



V(w, k) = J ± [T(w, k)f 



(11) 



where Y is given in Eq. (g). Note the absence in Fig. 2 
of gauge field lines connecting fermions on one plane to 
fermions on the other plane. This absence follows from 
the assumption of no coherent fermion hopping between 
the planes and is the reason why the interaction is en- 
hanced by the gauge field, rather than suppressed by it 
as is the in-plane interaction. 



' T T 



+ 



+. 



FIG 2. Ladder sum leading to pairing of fermions on 
adjacent planes. Solid lines denote fermions propaga- 
tors renormalized by the gauge field, dashed lines denote 
spin-spin interaction, J±, between the planes and solid 
triangles denote vertices renormalized by the gauge field. 



Because the interaction V(to, k) connects fermions on 
different planes it does not give rise to a fermion self 
energy in the leading order of perturbation theory. Di- 
agrams of higher order in Jj_ may be absorbed into the 
short-range interaction W between the fermions on each 
plane and its renormalization by gauge fields. We have 
previously shown that W is renormalized to zero by gauge 
fields, so these diagrams may be neglected. 

Therefore we may obtain the pairing effects of the in- 
teraction V from the ladder sum in Fig. 2. The analytic 
expression corresponding this diagram is, after integra- 
tion over pu 



E 



T J ± a 2 

47T Vp 



(12) 



A(p' ± +p±,e + uj)dp' ± 



VFP2 

Po^a 



3<r/2 



(w + e) 2 /3 w i/ 3 



The integral on the right hand side of this equation is 
infrared dominated and so may be evaluated by scaling 
(up to numerical factors). We also use the definition of 



wo 7 Eq- @, to eliminate the combination 
that T s is given by 

T s = (3uj (J ± a 2 g 2 )^ 



^ We find 

V F 



(13) 



where (3 is a numerical coefficient of the order of unity 
and g 2 is renormalized fermion-gauge field interaction 
constant. The value of g 2 has been estimated from the 
temperature dependence of the resistivityo, and u>o from 
the band structure. Roughly, l/(a 2 g 2 ) ~ 50 meV and 
ujo ~ 100 meV. Jj_ can be estimated from the cross re- 
laxation experiment to be Jj_ ~ 5 meV. In order for this 
T c to be of a reasonable order of magnitude, one must 
have 6(7 — 2 w 3 i.e. a « 5/6, because J± is small. This 
value of a would imply a 1/T%T ~ j 1 - 4 / 3 j a somewhat 
more rapid variation than is observed. 

Although the divergent interaction exists at all points 
on the Fermi line, the amplitude varies with position, 
9, because the gauge interaction energy scale luq varies 
withJ?. In fact, for the reasons discussed in Section VI of 
Ref.o we expect this variation to be substantial in high 
T c materials. Therefore we expect that as one lowers the 
temperature the gap first appears at the points 9 max at 
which ljq is maximal, so T s = /3u} (9 max )(J±a 2 g 2 ) 3 / 6cr ~ 2 . 
For T near T s the gap function A(9) will be very strongly 
peaked about 9 max ; indeed from Eq. ( |l2| ) we see that 
A($) decays away from 9-max as \9 9 max \ - ba . For T = 
the gap spreads over the whole Fermi surface but remains 
very anisotropic: 



A(0) w pu (e)(J x a 2 g 2 y 



(14) 



Finally, we note that because the gap function is so 
strongly peaked at particular points on the Fermi line, 
the energy is not very sensitive to the symmetry of the 
gap. In single-plane models, the pairing tends to be d- 
wave so that the two members of the Cooper pair can 
avoid each other. In the present model, the two members 
of the Cooper pair reside on different planes, so d-wave 
pairing becomes favorable only when tunneling between 
planes is includedo. 



V. GAP EQUATION: SPIN WAVES 

In this section we consider the implications for the for- 
mation of the spin gap of an alternative picture of the 
origin of the antiferromagnetic correlations. We suppose 
that there are weakly damped propagating spin waves n 



with dispersion ur 



! [(fc — Q) 2 ] + 5 2 and we ask how 



these lead to pairing. Here Q = (jr, n), c is the spin wave 
velocity and 5 is the spin wave gap. In the "quantum crit- 
ical" regime relevant for high T c superconductors, <5 = aT 
with awl. We assume that AG = \Q\ — 2pp > 0, so the 
low-energy spin waves lie outside of the particle- hole con- 
tinuum, and it is consistent to assume they are coupled 
to the fermions, but are undamped. The condition that 
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the low-energy spin-waves are outside the particle-hole 
continuum is (if c < vp) 



cAG > T 



(15) 



In high-T c materials we estimate from the Fermi line 
observed in Ref.E9 that AG_> O.lA -1 . In insulating 
La,2CuOi, c — 0.75 eV — JH; the previously discussed 
fits to T 2 /TxT in YBa 2 Cu 3 6 . 6 imply c = 0.35 eV-MB. 
Adopting the latter value we find that the condition for 
the validity of this assumption is T < 0.035 eV 400K. 
Thus the model may be adequate for discussion of spin 
gap phenomena at T ~ 150K, but the relevance to the 
spin dynamics at room temperature or above is question- 
able. 

The action describing the coupling of fermions to spin 
waves is 



Hf — H -\- H s 
Haw = ff y~^^, e, 

+ 51 TnU, 



(16) 



u),k 



2D(cj,k) 



-. -. j ^(1) ~(2) 



Here H is the fermion-gauge field action given in Eq. (]?]), 
n describes the undamped spin fluctuation and the spin 
fluctuation propagator D(uj,k) is 



D(u,k) 



u 2 + (cfc) 2 + S 2 



(17) 



and u and g are parameters with the dimension of energy. 
As in the previous section, the interplanc coupling leads 
to an effective pairing interaction V e f / given by 



V ef f(uj,k) 



(gu) 2 J ± 



(lu 2 + (ck) 2 + 5 2 ) 2 



(18) 



In order to account for the Cu relaxation rate data, we 
must assume that for T > T a the model is in the quan- 
tum critical regime in which 8 = aT. Note that the 
momentum k is measured from the commensurate anti- 
ferromagnetic wavevector Q = (it, 7r) and that V e ff is a 
very strongly peaked function of k. Thus, V e ff scatters 
a fermion of momentum p mostly to states with mo- 
mentum very near p + Q and if Eq. ( |l5| ) is satisfied 
and p is on the Fermi surface, these final states are far 
from the Fermi surface. This strong momentum depen- 
dence means that the gap equation decomposes into two 
equations, one expressing the gap at momentum near p, 
A e (p) , in terms of the Green function at momenta near 
p + Q and another one expressing the gap at momentum 
p + Q, A*(p), in terms of the Green function at momenta 
near p. Specifically, 



A ( ) = V Ty e//(^,g)A* + Jp + q) 

uT\e+^+~e p+q + [K + dv + qW 
A* (p) = V TV eff{u,q)A t+L0 (p + q) 



(19) 



Here C, p is the fermion energy, ( p = ( p +q; for circular 
Fermi line ( p = v F (\p\ — p F ). Because V e ff is a strongly 
peaked function of q, we may neglect the q in denomina- 
tor of ([l9]). The magnitude of the denominator is then 
controlled by the magnitude of Cp+Q- This depends on 
the position of p along the Fermi line. The minimal value 
of Cp+Q is v F AG and the maximal value is of the order of 
e F . Because J V e ff(uj,q)d 2 q is a singular function of lu, 
the frequency sum is dominated by to ~ max(T,d) and 
we expect that 5 ~ T, so by we may also neglect the 
frequency and A* dependence of the fermionic denom- 
inator in (19). The first equation ( |l9| ) then becomes a 
linear convolution equation; it can be combined with the 
second equation (19) giving an equation which contains 
only the gap function in the vicinity of the Fermi line 



Si 



W e ff(uJ,q)A e+UJ (p + q) 

^e+Jp + q) 



,.q — ^ 2/3 |e + ^/3 + C 2 +g 



(21) 



with the singular kernel 



W Bff (u,q) = Tj^Veffiv, k)V eff {Lo -rj,q-k) (22) 

ij.k 



In the gap equation ( |21| ) the assumption that p is on the 
Fermi surface means that the singularities in the fermion 
denominator must be also considered and compared to 
the singularities of the kernel W e ff- The fermion propa- 
gator depends sensitively only on the component of mo- 
mentum normal to the Fermi surface and, as we show 
below, A e (p) also changes smoothly along the Fermi sur- 
face. We may therefore integrate over the component 
along the Fermi line, q±, immediately, obtaining an equa- 
tion similar to (|2l|) but with q± = and a modified kernel 



1 



16c 2 ^ In 2 + (ckn) 2 + S 2 W 2 

TJ,fe|| 11 ' 

1 



[(b) - rj) 2 + (ck\\ - cq\\ ) 2 + 5 2 } 3 / 2 
Two cases then arise. If 

(^ /3 T 2 / 3 ,A) 



5 max 
- < - 

c Vp 



(23) 



(24) 



then the singularity in We// is dominant and we may 
perform the fen and gn integrals obtaining 



with 



T 



W eff {m)A e+ul (p) 



Q+q „ ^ 2/3 k + ^i 4/3 + e P + ^ + jp) 



(20) W eff (u) 



T(gu) 4 J 2 



16n 2 c 



[r] 2 + S 2 } [(uj - j]) 2 + S 2 } 



(25) 



(26) 
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However, if the condition rt24|) is not satisfied, then the 
singularity in the fermion denominator is dominant and 
we obtain 



A e (p) = 



T 



p+Q 



E 



W eff (u)A e+u (p) 



OJ, 



2/3, 



w| 4 /3 



with 



W eff {u) 



32ttc 3 vf 



E 



dk 



(27) 



[r; 2 + <5 2 + /c 2 ] 3 / 2 



[(w - r?) 2 + £ 2 + fc 2 ] 3 / 2 



(28) 



This system of equations may be solved to determine T s 
and A(T). At low T, or for not too small c we expect Eq. 
( pi| ) to be satisfied, so we shall consider in detail only the 
Eqs.( p5[]26| ). The other case leads to very similar results. 

To estimate T s we linearize Eq. (p5| ) and, because the 
sums are infrared dominated take only the contribution 
from the lowest Matsubara frequency. We find 



T s (p S ) 



1 3/5 



L47rV3 c 2 Cp()+QW V3 a2 j 



(29) 



Similarly, if the condition (|2J) is not satisfied we find 
from (p7]) and (E8h the slightly different formula 



T a (pe) 



v/3 



J±(gu) 



16^/3 c 3/2 l ,y 2 Cpe+QW l/6 a 5/2_ 



6/11 



(30) 



As found in the previous section, the onset of the spin 
gap is angle dependent, ft appears first at the angle, 
6*, for which Cp+Q has a minimum, and as T is low- 
ered spreads over the Fermi line. Because the interaction 
is more singular than in the case of the damped spin 
waves considered in the previous section, T s goes as a 
smaller power of J± . Roughly, A is large for 9 p such that 
T s (pe) > T, and drops rapidly for larger 9. We denote 
the interval in which A is large by (0* + 8 Q ,6* - 9 Q ). We 
note that 



(pe+Q — C* 



■ eo( 



(31) 



where eo is an energy scale of the order of cf- Thus, for 
T very near T s , 9q is given by 



T s -T 



1/2 



1/2 



For T much less than T s 

To 
T 



5/6 



1/2 



(32) 



(33) 



For \9 — 9* | > 6*o A(9) is induced by the value of A inside 
the interval (0* + 6 , 0* - 6 ). From Eqs. (|||6|) it can 
be seen that the kernel is so sharply peaked that 



A(0)~A(0 o )^- 
Finally, at T = one finds 

J±(gu) 2 



>) 4 



(34) 



A(9) 



47rc 2 Cp+Q^o /3 a 



3/5 



(35) 



In words, because the interaction is so strongly peaked, 
the value of the T — gap is controlled by the energy, 
(p+q, of the intermediate state of momentum p + Q. 



VI. CONCLUSION 



We have argued that any theory of the spin gap ob- 
served in Y BaiCu-sO^+x or Y Ba^Cu^O^ must be based 
on a pairing instability of a Fermi sea of chargeless 
fermions (i.e. a spin liquid) in the presence of antifer- 
romagnetic correlations. The large inelastic scattering 
found in the models of spin liquids leads to strong pair- 
breaking effects which can only be overcome by singular 
pairing interaction. The most plausible origin of this in- 
teraction is a between planes coupling enhanced by in- 
plane antiferromagnetic correlations. We considered two 
specific models: a spin liquid with 2pp over damped 
magnetic correlations induced by a gauge field interac- 
tion and a spin liquid coexisting with weakly damped 
antiferromagnetic spin waves. (A third possible model, 
namely a spin liquid with magnetic correlations induced 
by tuning-ainur fermion interaction, has been considered 
elsewhereE3o) The former model is clearly appropriate 
to optimally doped materials, where there is convincing 
photoemission evidence for a large Fermi line. We ex- 
pect by continuity that it is also appropriate for dopings 
somewhat below optimal. Weakly damped spin waves ex- 
ist in the insulating "parent compaund" materials, and 
it has been proposed that low doping induces a gap (but 
no damping) in the spin wave spectrum as well as a par- 
ticle hole continuum of spin excitations. This picture 
might be justified if there is a small Fermi line (i.e. "hole 
pockets" ) , but if there is a large (Luttinger) Fermi line 
the absence of spin wave damping is difficult to justify. 
It is controversial which picture applies to the most ex- 
tensively studied spin gap compounds Y BCL2CU3OQ1 or 
YBa,2Cu40s. We analysed both models, obtaining esti- 
mates for the onset temperature, T s , of the spin gap and 
the angular dependence of the gap function. 

In both models T s scales as a power of the between 
planes coupling J±, and the gap function is sharply 
peaked about particular regions of the Fermi line. In 
both models, pairing of spinons does not imply a true 
thermoduynamic transitioned. T s is a crossover temper- 
ature and superconductivuty sets in at a lower tempera- 
ture T c at which the charge carrying bosons condense. 
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In the overdamped case the physics is controlled by 
divergences in the fermion 2pp response function due 
to the singular gauge field interaction. These divergen- 
cies are characterized by an exponent a which depends 
only on the fermion spin degeneracy N and by an en- 
ergy scale loq defined in Eq. ([)]) which varies substan- 
tially as one moves along the Fermi line. We found that 

T s oc uj^ ax j]^ (uj 7 nax is the maximal value of lu (8) 
on the Fermi line which occurs at 9*). At T sa T s the 
gap is very sharply peaked at 9 = 9*. As T is lowered 
the region where gap is appreciable grows and the T = 

3 

gap A(0) oc Uo(9)Jx~ 2 ■ Evidently the result depends 
sensitively on the exponent a which has been estimated 
to be in the range 1 > a > 1/3. Consistency with NMR 
for T > T s requires that 1 > a > 2/3; for this range of 
a T s is proportional to J± to a power of order 1 and is 
of the correct order of magnitude, but the precise value 
depends on numerical coefficients which are not known. 

For T s > T > 0, the spin gap is appreciable over a 
part of the Fermi line, and suppresses the contribution 
of that part of the Fermi line to the uniform suscep- 
tibility and NMR relaxation rates. We see from Eqs. 
(PI [Til) that the gap is largest along the zone diagonal 
(where vp is largest) and smallest at the zone corners 
where vf becomes very small. All parts of the Fermi line 
make roughly equal contribution to the oxygen relaxation 
rate and uniform susceptibility (although the logarith- 
mic divergence associated with the van Hove singularity 
may emphasize the corners to some extent), so we may 
roughly estimate the suppression of these quantities from 
the fraction of the Fermi line which is gapped. The cop- 
per relaxation rate is more complicated. It is dominated 
by the 2pp fluctuations which lead to 



Thus the contribution to the Cu relaxation rate is 
largest where u>o is largest (i.e. along the zone diago- 
nal) so one would expect that in this model the forma- 
tion of the spin gap would suppress the Cu relaxation 
rate more strongly than the oxygen rate. However, the 
contributions of the ungapped portions of the Fermi line 
continue to grow as T is lowered, so the maximum in the 
Cu 1/T\T relaxation rate occurs at a T < T s determined 
by the interplay between these two effects. 

We now consider the underdamped case. The basic as- 
sumption is that there are two distinct types of spin ex- 
citations: propagating antifcrromagnetic spin waves and 
particle-hole continuum of spinon excitations. This pic- 
ture has been derivedo from a microscopic Hamiltonian 
using the assumption that the doping is so low that long 
range magnetic order is present, and it is plausible that 
it may apply to lightly doped high T c materials which 
lack long range order if these materials do not have large 
(Luttinger) Fermi line but have instead "hole pockets" . 
An advantage of this picture is that large q properties 



are dominated by spin waves, which explain in a natural 
way the strong T dependence of the Cu 1/T\T and I/T2 
rates observed experimentally. The disadvantage of this 
picture is that the same spin waves would give a factor of 
six too large contribution to dx/dT at high temperatures 
so to account for the observed dx/dT one must assume 
that the particle-hole continuum leads to a large negative 
contribution to dx/dT which almost precisely cancels the 
spin wave contribution. 

We studied the pairing of fermions in the presence of 
the undamped spin waves. The pairing interaction is 
very strongly peaked at the antifcrromagnetic wavevector 
Q = (tt, 7t), so the dominant process in the gap equation 
is a virtual scattering from a state p on a Fermi line to 
one at p + Q away from the Fermi line. We found that 
the resulting pairing interaction is very strong, the on- 

3 /5 

set temperature T s oc J ^ . As in the overdamped case, 
the gap has a strong angular dependence, it appears first 
at a particular point on the Fermi line and spreads over 
it as T is decreased. However in the underdamped case 
the angular dependence of the gap is controlled by the 
energy of the intermediate state, Cp+Q- The pairing af- 
fects the particle-hole contribution to physical response 
functions but does not directly affect the spin wave con- 
tribution. Thus, the copper relaxation rate is not signifi- 
cantly affected by the pairing, but the oxygen relaxation 
rate and uniform susceptibility are. The detailed tem- 
perature dependence is determined by the way in which 
the gap spreads over the Fermi line as the temperature is 
decreased, and this depends sensitively on the shape of 
the Fermi line and, in particular, on its curvature. 

This picture provides a natural explanation of the 
magnetic properties of lightly doped YBaCuO materi- 
als (Y Ba2Cu?,0§+ x with 0.45 < x < 0.6) in which the 
T s inferred- from the uniform susceptibility is of the or- 
der 300-ftTE3 while strong antiferromagnetic fluctuations 
with no evident gap have been observed in neutron 
scatteringpH. However, in this picture the coincidence 
of T s with the maximum of the Cu 1JT\T rate observed 
in Y BaiCu-iO%,i and Y BaiCu^O-g, is an accident. 

Thus, we believe that the overdamped picture is more 
appropriate to materials with doping greater than one 
in YBCI2CU3O6.J, while the underdamped picture may 
be more appropriate for more lightly doped materials. 
This general point of view has been advocated by igwy 
authorstS and has been emphasized recently in Ref.BE£l. 
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